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We compute low-lying eigenmodes of the gauge covariant Laplace operator on the lattice at finite temperature. 
For classical configurations we show how the lowest mode localizes the monopole constituents inside calorons and 
that it hops upon changing the boundary conditions. The latter effect we observe for thermalized backgrounds, 
too, analogously to what is known for fermion zero modes. 

We propose a new filter for equilibrium configurations which provides link variables as a truncated sum involving 
the Laplacian modes. This method not only reproduces classical structures, but also preserves the confining 
potential, even when only a few modes are used. 
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1. Introduction 

As an analyzing tool we study low-lying eigen- 
modes of the gauge covariant Laplace operator 

4 

-A(/)„(x) = Xn(t>n{x) (1) 

with Ufj_{x) a given lattice configuration in the 
fundamental representation of SU(2). Compared 
to fermionic (near) zero modes, Laplacian modes 
do not suffer from chirality problems nor dou- 
bler problems. They seem not related directly to 
topology (like zero modes are, due to index the- 
orems), but it has been observed that they are 
sensitive to the location of instantons |1I2| . The 
scaling properties of their localization have been 
investigated recently 

Here we will present two ideas concerning 
Laplacian modes @]. The first one is to study 
what can be learned from the profile of the mod- 
ulus of the lowest mode (the 'ground state proba- 
bility density') \(/)o{x)\'^, in particular with chang- 
ing boundary conditions. The second idea is to 
introduce a novel low-pass filter. It concerns the 
reconstruction of the gauge background from rela- 
tively few Laplacian modes. The purpose of both 



methods is to identify the underlying infrared de- 
grees of freedom in the gauge field that are re- 
sponsible for features like confinement. 

2. Profile of the lowest mode 

2.1. Caloron backgrounds 

As a testing ground for the Laplacian modes 
we first investigate a caloron of maximally non- 
trivial holonomy with its monopole constituents 
|5I6| put on a 32^ • 4 lattice^. Fig. [T] (top) shows 
the action density along the line connecting these 
monopoles, which are clearly visible as two self- 
dual (and almost perfectly static) lumps. As 
another signal the Polyakov loop goes through I2 
and — 12 - which amounts to a local symmetry 
restoration - at the monopole cores. 

The lowest-lying mode of the Laplacian in this 
background is shown in Fig. (bottom) IJ. One 
can see that the presence of monopoles is reflected 
by a maximum resp. a minimum in the profile 
of the lowest mode (although slightly shifted). 
In addition, we have depicted the lowest-lying 
Laplacian mode with antiperiodic boundary con- 
ditions in Euclidean time. In that mode the 
role of the monopoles is interchanged, which can 
be understood from a symmetry of the caloron. 
Hence the lowest-lying Laplacian mode 'hops' as 
the result of changing the boundary conditions. It 
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^This is done by calculating links from the continuum 
gauge field, followed by a few steps of cooling. 
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Figure 2. The kinetic term for the caloron of 
Fig. □ with periodic (boldface) and intermediate 
(dashed, muhiphed by 8) boundary conditions. 
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Figure 1. Top: action and topological density 
(not distinguishable, both multiplied by a factor 
500) and the traced Polyakov loop shown along 
the line connecting the constituents in a large 
caloron on a 32^ • 4 lattice. Bottom: profile of 
the lowest Laplacian mode in that background 
with periodic (boldface), antiperiodic (thin) and 
intermediate (dashed) boundary conditions, re- 
spectively. 



behaves similar to the caloron fermion zero mode 
[7] in the context of which complex boundary con- 
ditions were introduced first. Also for the Lapla- 
cian modes we now allow for a general phase in 
the boundary condition 

(j){x4 + Ni) = exp(27rO0(x4) . (2) 

The angle ^ can be restricted to C G [0, 1/2] since 
the remaining ^'s can be mapped to this inter- 
val by charge conjugation. For the intermedi- 
ate case C = 1/4, i.e. halfway between periodic 
and antiperiodic boundary conditions, we find a 
valley in the mode profile extending between the 
monopoles (dashed curve in Fig. (bottom)). 



The Laplacian eigenmodes are stationary 
points of the functional T[Uf^: <j)] — J^x ^(^) with 

tix)^Y.p;:''ix)q^''{x+f,)-rix)\' (3) 

the lattice analogue of the square of the covariant 
derivative. This kinetic term (in the notion of 
Higgs models) represents the local distribution of 
the remaining 'energy' when the mode (j) is best 
adapted to the background U^. Here we view 
the quantity t{x) as another observable that is 
sensitive to structures in the lattice configuration 
under consideration. 

Fig. 121 shows the kinetic term for the caloron 
and its lowest mode with different boundary con- 
ditions. With the periodic mode the kinetic term 
reveals the monopole by a volcano-like structure, 
i.e. a maximum with a dip in its center. Like 
the mode itself the kinetic term hops to the com- 
plementary monopole when changing boundary 
conditions to antiperiodic (not shown). For the 
lowest mode with intermediate boundary condi- 
tions both monopoles are detected by maxima in 
this observable (dashed curve in Fig. (SJ. 

2.2. Thermalized backgrounds 

In this subsection we explore the Laplacian 
modes as an analyzing tool on thermalized gauge 
field configurations 0]. Fig. |31 shows the lowest- 
lying mode in a generic background obtained on 
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Figure 3. First three rows: the lowest Laplacian mode in a thermahzed background with C = 0.02, 0.355, 
0.48 from left to right and in lattice planes {x2,X4) = (12, 1), (4, 4), (6, 3) from first to third row (vertical 
scale 0.05). The last three rows show the corresponding profiles of the kinetic term (vertical scale 0.0025). 



4 



a 16'^ ■ 4 lattice at f3 — 2.2 (which amounts to 
T = 0.75 Tc) confirming the expectation that 
these modes are free of UV fluctuations. The 
'evolution' of the Laplacian mode with the an- 
gle C is shown horizontally in Fig. |3 while the 
first three rows represent different but fixed lat- 
tice planes, which contain the global maximum of 
the mode. In the example we find three lattice lo- 
cations which become the global maximum within 
some ^-interval. Therefore, the lowest Laplacian 
mode in a thermalized background hops, too. 

At the hopping points in C the value of the max- 
imum^ as well as the inverse participation ratio 
(taking values between IPR=4 and IPR= 11) de- 
crease^ and the mode rearranges itself. For some 
boundary conditions the lowest Laplacian mode 
can be characterized as a global structure @] . As 
the three last rows of Fig. |3| show, the kinetic 
term follows the modulus of the mode, but is less 
smooth. 

Analyzing 50 independent configurations we 
found that the lowest mode hops up to four times 
as a function of C. Apart from the Laplacian 
mode being wave-like, the hopping phenomenon 
resembles the behavior of the fermion zero mode 
in thermalized configurations [S]. However, in 
both cases it is not obvious which gluonic feature, 
e.g. in the topological charge, distinguishes the 
preferred locations. Interestingly, we have found 
a correlation of the maximum of the periodic and 
antiperiodic Laplacian mode to positive and nega- 
tive Polyakov loops, respectively, see Fig.^ This 
is the same tendency as for calorons (cf. Fig. ^ 
which could be understood as an argument in fa- 
vor of calorons 'underlying' the quantum config- 
urations. Alternatively, the Polyakov loop could 
play a role defining pinning centers for the Lapla- 
cian mode in the spirit of Anderson localization 
in a random potential 

The lowest Laplacian mode in the adjoint rep- 
resentation we have found to have minima at the 
caloron monopoles (see also llllOp that extend 
to two-dimensional sheets for antiperiodic bound- 
ary conditions In thermalized backgrounds. 



periodic 
antiperiodic 



^The global minimum is not stable enough to employ it 
for analyzing purposes. 

''The C- values of Fig. El are chosen such that the IPR is 
maximal within the corresponding (^-interval. 
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Figure 4. Scatterplots of the Polyakov loop (hor- 
izontal) at lattice sites where the modulus of the 
lowest Laplacian mode (vertical) is large: stars 
represent the periodic mode, circles the antiperi- 
odic one. 



the maxima of the adjoint modes are correlated 
to the fundamental ones with the same bound- 
ary condition, but are stronger localized than the 
latter. 

3. A Fourier-like filter 

Laplacian ('harmonic') eigenmodes can be used 
to define a new Fourier-like filter which we present 
now 4?. We were inspired by the representation 
of the field strength in terms of fcrmionic modes 
in but here we aim to reconstruct directly 
the link variables on which any observable can be 
measured. To this end we combine the definition 
of the lattice Laplacian, Eq. with a spectral 
decomposition and at y — x + (1 immediately ob- 
tain 

AT 

JV ^ N1N2N3N4 ■ 2 . (4) 

The idea is to truncate the sum on the right hand 
side of this equation at a small mode number N 
compared to their total number Af. 

The question arises, how to relate such an ex- 
pression to a unitary link variable. Here the 



charge conjugation helps, since it guarantees 
that every eigenvalue is two-fold degenerate with 
eigenfunctions related as (/>'"' {x) = e"''' (/)*'' (x) . The 
corresponding bilinear expressions in Eq. add 
up to an element of SU{2) up to a factor. Such 
a situation is familiar from the staple average in 
cooling or smearing. We divide by the square 
root of this factor, which is the determinant and 
positive for all practical purposes, and obtain the 
final filter formula 

N 

e-e(^)0'(^ + A)e"'])/yd^^rO- (5) 

The quality of the filter is controlled by N, 
where N= Af reproduces the original configura- 
tion exactly. In the other extreme case, N = l, 
the filtered links U^{x) can be shown to be 
pure gauge. From the behavior of 4>{x) under 
gauge transformations it follows immediately that 
Ufj,{x) transforms as a link (and no gauge fixing 
is involved in the filter). 

The filtering can be performed with Laplacian 
modes of any boundary condition. Then Eq. jSJl 
is just the average over opposite boundary con- 
ditions C and — C. This additional parameter C 
completely fixes the Polyakov loop at iV = 1 to 
tr'P(x)/2 = cos(27rC), while for nontrivial cases 
N > 2 the Polyakov loop is observed to fluctu- 
ate around this value. It follows that in order 
to optimize the filter for our circumstances (con- 
fined phase, nontrivial holonomy calorons) the 
best choice is the intermediate case ( = 1/4. 

In order to test the filter, we again start with 
a caloron (this time over 16^ • 4). The action and 
topological density as well as the Polyakov loop 
plotted in Fig. [S] show, that the filter with the 
number of modes as low as A'^ = 4 starts to repro- 
duce the classical structures qualitatively, while 
for N — 150 the agreement is almost perfect. 

During the filtering we have recorded the de- 
terminant used in Eq. © to project to an SU{2)- 
element. With increasing number N of modes it 
increases, too (i.e. becomes closer to 1). More- 
over, its variation over the lattice decreases, such 
that all links appear 'equally well filtered'. 
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Figure 5. Action and topological density and 
Polyakov loop of a large caloron on 16"^ • 4, fil- 
tered with boundary condition C = 1/4. From 
top to bottom: = 4 modes (densities multi- 
pHed by a factor of 100), N — 150 modes (factor 
400) and the original configuration (factor 400). 



Now we turn to the application of the filter 
to thermalized configurations (at finite temper- 
ature), first by the example of Sect. 12.21 

The first question we want to address is the 
correlation of the filtered links to the original 
ones. Fig. El shows scatterplots including all 
SU (2)-components of the links in all /i-directions 
and over all lattice locations. To compare link 
variables themselves makes sense, since the orig- 
inal and the filtered configuration come in the 
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Figure 6. Correlation of link variable components txU^{x)/2 and trU^{x)Ti/2i [ri the Pauli matrices) 
between the original configuration and the filtered one with N = 2 (left) and N = 100 (middle) and 
within filtered links with = 2 vs. iV = 4 (right). In order not to overload the plots, only 2000 link 
variable components have been plotted. 



same gauge (in other words the application of 
the filter commutes with gauge transformations). 
Fig. El (left) shows that even in the first nontrivial 
case N — 2 the filtered links are correlated to 
the original ones. With increasing number of 
involved modes this correlation becomes stronger 
as expected: in the scatterplot for N — 100 
(Fig. El (middle)) the deviation from the diagonal 
decreased compared to iV = 2. One can make 
this statement quantitative by computing the 
average square distance from the diagonal, which 
indeed is smaller for bigger N. Actually, it also 
shows that the time-like links are less correlated 
than the space-like ones, especially when using 
the 'inappropriate' periodic boundary conditions. 
Finally, Fig. (right) shows that the filter with 
comparable number of modes produces strongly 
correlated links. 

What is even more remarkable is that the filter 
method preserves the string tension. We plot the 
logarithm of the Polyakov loop correlator mea- 
sured on 50 configurations at T = 0.75 Tc in Fig. 
13 (left) , both for the original one and the filtered 
ones with N = 2, 10, 100. It reveals a clearly lin- 
ear behavior with a slope reproducing the original 
one within 15%. The minimal number of modes, 



N — 2, suffices for this behavior. We conclude 
that the confining properties of lattice gauge the- 
ory are captured by the lowest Laplacian modes. 

This property actually does not depend much 
on the ordinal number of the modes, see Fig. [71 
(right) . First of all, one can include the N highest 
Laplacian modes. Due to the symmetry [31 

0.V-n = (-1)^^ ^"'/'n , Aa/'_„ = 16 - A„ , (6) 

including these highest modes results in the same 
bilinear contributions in Eq. (0) with the only 
modification A„ 16 — 2A„. The Polyakov loop 
correlator computed from the 50 configurations 
filtered in this way agrees very well with the one 
discussed so far (and is therefore not shown in 
Fig. I?!). The second possibility is to shift the 
mode number from 1 till 4 to say 5 till 8 or 51 
till 54. The corresponding curves shown in Fig. 
H (right) reveal the same linear behavior up to a 
vertical shift. The Polyakov loop correlator is ro- 
bust also against a third option, namely replacing 
the eigenvalues A„ in Eq. © by a constant, which 
again produces identical curves (not shown). 

After filtering the Polyakov loop correlator 
has no sign of a Coulomb regime, since the fil- 
ter has washed out short range fluctuations. At 
zero distance the values {{trV(x))'^) are smaller 
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Figure 7. Polyakov loop correlators of 50 configurations at finite temperature, filtered with boundary 
condition C = 1/4 and plotted on an inverse logarithmic scale over the distance. Left: the original 
configurations vs. the filtered ones with N — 2, 10, 100 modes. Right: 4 modes with various shifts in the 
ordinal number of the modes. 



than the original one. This reflects the already 
mentioned effect that the filtered links for small 
TV are mostly traceless (governed by the choice 
( = 1/4). The site distribution of the Polyakov 
loop is a narrow one around trP = for N = 2, 
but for N = 100 is compatible with the original 
distribution which in turn is described by the 
Haar measure. Fig. |S1 shows locally (for the ex- 
ample configuration and in a fixed lattice plane) 
that with growing N more and more fluctuations 
appear in the Polyakov loop. Studies are under 
way to clarify whether the new filter method 
preserves other infrared features, like the spatial 
Wilson loops, and how it behaves at zero tem- 
perature. 

Since the filter preserves the string tension from 
Polyakov loop correlators and is not biased to 
classical solutions nor to particular degrees of 
freedom in the gauge field (like monopoles or vor- 
tices), it is interesting to have a closer look at the 
vacuum structures that emerge when the filter is 
applied to generic configurations. 

Concerning the action density, these objects 
are found to be isolated peaks (which are non- 
static and not necessarily (anti)selfdual). This 
phenomenon seems counterintuitive as the filter 



uses the lowest-lying Laplacian modes which are 
smooth. However, too small action density lumps 
have also been observed when filtering smooth 
calorons at low N (cf. Fig. |5l (top)). This phe- 
nomenon presumably comes from the nonlinear- 
ity in the filter formalism (Eq. Q), in particular 
the determinant used to scale up the mode con- 
tributions to an 5L/(2)-element. Interestingly, we 
have found a correlation of the peak structure of 
the filtered action density to that emerging after 
cooling or smearing in an early stage. 

We stress that, in contrast to e.g. cooling, the 
filter does not involve the plaquette action and 
thus is not made to reduce it in the first place. 
Upon filtering the links, the total action does be- 
come smaller, but only down to some amount 
(in our example down to 66 instanton units at 
= 2), which might well be the minimal con- 
tent necessary to reproduce the infrared physics. 
More work has to be done to better understand 
the spiky structures induced by the filter. 
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Figure 8. Polyakov loop 'evolution' with increas- 
ing number of modes for the thermalized exam- 
ple in a fixed lattice plane {x2 = 6). From top to 
bottom: N = 2, 10, 100 and the original configu- 
ration. The vertical axes are from —1 to 1. 
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